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Opposite-sign lepton pairs (ee,e/i and /i/i) from open charm decay are proposed as a 
measure of nuclear shadowing effects. Via an approximate scaling the ratio of the dilepton 
<^ spectra from p + A to those from pp reflects the shadowing function well. We show that 

■ the required measurements are feasible at the Relativistic Heavy ion Collider (RHIC) by 
considering the backgrounds according to the proposed PHENIX detector geometry. 
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1.1 Dileptons as Signals 



1. INTRODUCTION 

> 
O 

The spectrum of dileptons produced in heavy ion collisions at high energies has been 
proposed in the past to provide information about the dynamical evolution of quark- 
gluon plasmas. Unlike hadronic probes, that spectrum is sensitive to the earliest moments 
in the evolution, when the energy densities are an order of magnitude above the QCD 
confinement scale (1 GeV/fm 3 ). For example, dileptons from Drell-Yan process probe the 
quark degrees of freedom, and thermal dileptons tell us about the characteristics such as 
temperature and phase transition in the dense plasma[I|]. However, the small production 
cross sections for thermal dileptons with invariant mass above 1 GeV together with the 
large combinatorial background from decaying hadrons necessitate elaborate procedures 
to uncover the signal from the noise. The PHENIX detector 0, now under construction 
at RHIC, is designed to measure ee, e/i, and /i/i pairs to carry out this task. One of the 
important sources of background in the few GeV mass range arises from semileptonic 
decay of charmed hadrons (DD). As shown recently in ref. the expected thermal 

signals in that mass range may only reach 10% of the number of pairs from open charm 
decay. Therefore, special kinematical cuts and precise e/i measurements will be needed 
to uncover the thermal signals. 

There have been several attempts || to take advantage of the large open charm back- 
ground as a probe of the evolving gluon density. Most mid-rapidity charm pairs are 
produced via gluon fusion. Therefore, the dileptons from charm decay carry information 
about the distribution of primordial gluons before hadronization. In fact, the inside- 
outside cascade nature of such reactions greatly suppresses all sources of charm produc- 
tion except the initial perturbative QCD source. Therefore, the open charm background 
is dominated by the initial gluon fusion (gg — > cc) rate, and it depends sensitively on the 
nuclear gluon structure function, (/^(x, Q 2 ). 
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1.2 Nuclear Shadowing Effects 

The quantity of fundamental interest is the gluon shadowing function 



Rb/a{x, Q 2 ) = 9a(x, Q 2 )/Ag N (x, Q 2 ). (1) 

Shadowing of the quark nuclear structure functions, qA^x, Q 2 ), is well established from 
deep inelastic iA — > £X reactions^. For heavy nuclei, the quark structure functions are 
shadowed by a factor, R q /A(x <C 0.1, Q 2 ) ~ 1.1 — 0.1 A 1 / 3 , nearly independent of Q 2 . 
Perturbative QCD analysis]!, || predicts that the gluon structure will also be shadowed 
at x < 0.01 due to gluon recombination processes. A systematic measurement of gluon 
shadowing is of fundamental interest in understanding the parton structure of nuclei. The 
gluon nuclear structure is also of central importance in the field of nuclear collision since 
it controls the rate of mini-jet production that determines the total entropy produced at 
RHIC and higher energies. Recently there appeared an attempt [|K|] to extract nuclear 



gluon shadowing from the high statistics FA data on S n and C. Similar to quarks, gluons 
showed shadowing at small x and anti-shadowing at medium x. 

1.3 Our Purpose 

The point of our studyjil] is that the A dependence of continuum dilepton pairs in 



the few GeV mass region provides a novel probe of the unknown gluon structure. We also 
show that the required measurements of ee, e/i, and /z/i pair yields in p + A — > £ + £~X 
at y/s = 200 AGeV are not only experimentally feasible at RHIC but also that the open 
charm signal can be easily extracted via the proposed PHENIX detector. We focus here 
on p+A collision rather than A+A to minimize the combinatorial 7r, K decay backgrounds 
and other final state interaction effects. 

A key advantage of the continuum dilepton pairs from open charm decay over those 
from J/ip decay is that it is possible to test the applicability of the underlying QCD 
dynamics at a given fixed y/s by checking for a particular scaling property discussed be- 
low. In quarkonium production the mass is fixed and the required scaling can be checked 
only by varying the beam energy. In p + A —>■ J/ip production the required scaling was 
unfortunately found to be violated in the energy range 20 < y/s < 40 AGeV[[L2|] thus pre- 
cluding a determination of Ra- Possible explanations for the breakdown of QCD scaling 
for J/ip production and its anomalous negative xp behavior [12, 13] include interaction 



of next-to-leading order quarkonium Fock state with nuclear matter |14j], nuclear and co- 



mover J/ip dissociation|l5[, and parton energy loss mechanisms [16[|. It remains an open 



question whether the required scaling will set in at RHIC energies 60 < ^fs < 200 AGeV. 

2. DILEPTONS FROM OPEN CHARM 

2.1 Shadowing in Dilepton Spectra 

To demonstrate the sensitivity of open charm dileptons to gluon shadowing R 9 /a(x, Q 2 ), 
we calculate the pair spectrum from open charm decay in p + Au reactions at 200 GeV/A 
using for illustration two different gluon shadowing functions. The first shadowing sce- 
nario assumes that gluon shadowing is identical to the measured quark shadowing, i.e. 
R g /A(x,Q 2 ) = R q /A(x,Q 2 = 4 GeV 2 ), and is thus essentially Q 2 independent. This is 
the default assumption incorporated in the HIJING Monte Carlo model |17|. The sec- 
ond is taken from Eskola[|l^], where R 9 /a(x,Q 2 ) is computed using the modified GLAP 



evolution || starting from an assumed R 9 /a consistent with the measured quark shadowing 
function at Q 2 = 4 GeV 2 . This second shadowing function depends on scale and differs 
for q and g. The two shadowing functions are shown in Fig. 1. 

The initial charm pair distribution from B + A nuclear collisions at impact parameter 
b is computed as: 

dN BA (b) _ T ,r»» sm , t =. ^v- t, „2^,da, 



dp 2 ± dy 1 dy 2 



-n - r a )Y,x b T b/B (x b , Q 2 , r b )x a T a/A (x a , Q 2 ,r a ) — 4- (2) 

J b.a dt 



where T a /A(x, Q 2 , r) = TA(r)f a /N(x, Q 2 )Ro,/a(x, Q 2 , t) is the nuclear parton density func- 
tion in terms of the known nucleon parton structure functions f a /^(x 1 Q 2 ) 1 the nuclear 
thickness function T^(f) = / dznA^sl z 2 + f 2 ), and the unknown impact parameter de- 
pendent shadowing function R a / A (x, Q 2 ,r). The conventional kinematic variables x a>b are 
the incoming parton light cone momentum fractions, t = — (p b — Pi) 2 , and the final c, c 
have rapidities y 'i,y 2 and transverse momenta pi± = —p 2 L = P±- Next-to-leading order 
corrections to the lowest order parton cross sections da ab /dt are approximately taken into 
account by a constant K factor. 

We use the recent MRSA[^] parton structure functions. From a fit to low energy 
open charm production data in pp collisions, we fix m c = 1.4GeV, K = 3, Q 2 = s/2 for 
gg — > cc and s for qq — > cc. This choice of parameters leads to a charm pair cross section 
of 340/X& for pp collisions at RHIC. For impact parameter averaged collisions, the integral 
over the transverse vector in eq.@ leads to a factor: BAR b /B{x b ,Q 2 )R a /A{.x a) Q 2 )/af r A , 
where af^ is the inelastic B + A cross section and R a /A( x a, Q 2 ) is the impact parameter 
averaged shadowing function. 

2.2 Hadronization Scheme for Charm 

In order to compute the DD pair distribution, a hadronization scheme for c —>■ D must 
be adopted. In e + e~ — > cc — > DDX hadronization can be modelled via string fragmenta- 
tion or fitted for example via the Peterson fragmentation function (see |[20|| ). The final D 
carries typically only a fraction ~ 0.7 of the original c momentum. However, in pp colli- 
sions charm hadronization is complicated by the high density of partons produced during 
beam jet fragmentation. In this system recombination or coalescence of the heavy quark 
with comoving partons provides another mechanism which in fact seems to be dominant at 
least at present energies. The inclusive pp — > DX data is best reproduced [p0|, pTj] with the 



delta function fragmentation D(z) = 5(1 — z) in all observed Xf and moderate px regions. 
This observation can be understood in terms of a coalescence model if the coalescence 
radius is P c ~ 400MeV. We assume that hard fragmentation continues to be the domi- 
nant mechanism at RHIC energies, and hence no additional A dependent effects due to 
hadronization arise. However, this assumption must be tested experimentally. Our main 
dynamical assumption therefore is that the pp — > D transverse momentum distributions 
can be accurately reproduced from the QCD level rates assuming hard fragmentation as 
at present energies ||21|| . At RHIC this can be checked either via single inclusive leptons|| 
or directly via Kit in STAR experiment. 

With this assumption, the impact parameter averaged DD pair distribution in p + A 
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Figure 1: The dilepton dN/dMdy ratio curves (shadowed over unshadowed) at different 
masses obtained via Monte Carlo calculation are compared with the shadowing curves. 
The upper solid curve is EskolafT8|j gluon shadowing for Q 2 = 10 GeV 2 , and the lower 
solid curve is from HIJING [O- The ratio curves are plotted as a function of the scaling 



variable log 10 xa from eq.(|). The scaled dilepton ratios reflect closely the input shadowing 
functions. 



is given by 

, 2~T A = f 3 E J Z X bfb/N(x b , Q 2 )x a f a/N (x a , Q 2 )R a /A(x a , Q 2 )^ (3) 

dp±dy 3 dy A af* E 1 E 2 ^ ' dt 

From this distribution of charmed mesons, we computed the dilepton spectrum via Monte 
Carlo using JETSET7.3 to decay the DD. 

2.3 Scaling 

Since gluon fusion dominates, the light cone fraction xa that one of the gluon carries 
from nucleus A is related to the dilepton observables M and y. So there is an approximate 
scaling: 

\nx A = -y c5 + ln(M c5 /v^) w -y + ln[(/3M + M )/^/7s] (4) 

In the above we have related the invariant mass and the pair rapidity from cc to the 
dilepton from that cc decay. While the lepton pair mass M and rapidity y fluctuate 



around a mean value for any fixed M c5 and y c5 , on the average (y) ~ (y C c), and (M cg ) 
can be well approximated by a linear relation (M c5 ) w /3M + M over the lepton pair 
mass range 1 < M < 4GeV. For the delta function fragmentation case, (3 ~ 1.5 and 
Mq « 3.0GeV are determined by the D-meson decay kinematics. 

Therefore the ratio of the dilepton dN/ dMdy spectra in pA to scaled pp for different 
pair masses is expected to scale approximately as 



where v = Acrf^/af^ ~ A 1 ^ . In order to test gluon dominance and the accuracy of the 
above approximate scaling we compare in Fig. 1 the gluon shadowing function to the above 
dilepton ratio from the Monte Carlo calculation. The solid curves are the input gluon 
shadowings as a function of the Bjorken x. The other six curves are ratios of dilepton 
dN/dMdy spectra of shadowed p + Au over those from unshadowed p + Au as given by 
eq.(||) as a function of the scaling variable xa for three different dilepton masses. E.g., in 
Fig. 1 we first plot all the ratio curves in terms of reversed pair rapidity —y, then shift 
the ratio curves at M — 1, 2, 4 GeV to the left by 3.79, 3.51, 3.10 respectively. 

Overall, the scaled ratio of lepton pair spectra approximates the shadowing function 
remarkably well. We conclude that this ratio can therefore serve to map out gluon shad- 
owing in nuclei. Note that in the case of Eskola shadowing, even the Q 2 dependence of 
the shadowing function is visible through the rise of the ratio curves with M in the small 
x region. 

2.4 Signal to Background in Dilepton Spectra 

It is important to estimate the dilepton background to see if the proposed signal is 
experimentally feasible. Therefore we also calculated the signal and background dilep- 
tons for the PHENIX detector]^] taking into account the detector geometry and specific 
kinematical cuts. The electron background is mainly due to tt° Dalitz and photon con- 
versions. The background muon arises from random decays of charged pions and kaons. 
The electrons from 7r° Dalitz decay can be suppressed by small angle cut on dielectrons, 
and the background muons are mainly suppressed by reducing the free decay volume. For 
the detector geometry, the electron arms cover pseudo-rapidity range —0.35 < r] e < 0.35, 
and azimuthal angle range ±(22.5°, 112.5°). The muon arm covers pseudo-rapidity range 
1.15 < r/ M < 2.44 and almost full azimuthal angle. For the kinematical cut, we take 
E e > 1 GeV, > 2 GeV, and require that the relative azimuthal angle of the lepton 
pair <p£+e- > 90° in order to improve the signal-to-background ratio. 

For this background study, we use central p+Au collision with HIJING shadowing. We 
calculate the ee, e/i, and /i/i signal and backgrounds coming into the detector, where the 
background electrons and muons are generated from HIJING Monte Carlo calculation. 
We find (see Fig. 2) that the signal-to-background ratio for ee is very large, thus the 
dielectron signal from open charm decay is the easiest to extract. That ratio falls to 
about 2.5 for e/i, and about 1/4 for /i/i. The pair rapidity distributions of the signal 
and backgrounds for central p + Au collision are shown in Fig. 2. Due to the detector 
geometry, ee, e/i, and /i/i spectra cover pair rapidity regions centered at about 0,1, and 2. 



K?£ (M,y = - \nx A + In [(/3M + Af )A/i] ) 





Figure 2: Charm signal and backgrounds coming into the detector are plotted as a 
function of the lepton pair rapidity y. The solid curves are the signals from open Charm 
decay — labeled CC signal. The label C refers to the lepton coming from open Charm 
decay, label D refers to the electron coming from Dalitz and photon conversion, and label 
R refers to the muon coming from Random decay of pions and kaons. In (a), the dielectron 
backgrounds are the dashed curve DC and the dot-dashed curve DD; and both are scaled 
up by a factor of 100. In (b), the opposite-sign e/x backgrounds are the dotted curve CR, 
the dashed curve DR, and the dot-dashed curve DC. In (c), the dimuon backgrounds are 
the dashed curve RR and the dot-dashed curve RC. 

With the addition of the second muon arm, one can measure pair rapidity regions around 
— 1 and —2. Like-sign subtraction should significantly reduce the noise, especially in the 
fifi channel. We conclude that the proposed measurement is feasible. 

2.5 Uncertainty from Energy Loss and Multiple Scatterings 

The Cronin effect and energy loss are estimated to lead to distortions up to 10% 
at RHIC energies. From studies of the nuclear dependence of transverse momentum of 
J /if) production, the typical increment, 5p\(A), due to multiple collisions is limited to 
~ 0.34GeV 2 even for the heaviest nuclei^]. This momentum spread is shared between 
the c and c, and distributed approximately as g (Sp ± ) = e^^/nA 2 , A 2 ~ 0.17GeV 2 . 
For energy loss, we assume that the charm quark loses an energy 5E in the lab frame, and 
the charm quark p± is reduced by (1— e), where e = 5E/[m± smh(y+yo)]. Combining these 

two effects, we may write for the charm quark that the final p± = (p l ™+5p±)(l— e). There- 
fore the D-meson spectrum F(p±) = d 2 N/dm\ becomes F'(p±) = J F(p t ™)g(5p±)d5p±. 
Taking F(p±) oc e~ apA -, a ~ 1.3 GeV" 1 from HIJING, and expanding the convolution to 
lowest order, the relative change of the D-meson spectrum is given by 



F'{p±)/F(px) - 1 - ap ± e + a 2 A 2 /4 



(6) 



This is also the relative change of the cc pair spectrum when M cg = 2m± in case of y\ = 1/2, 
p±i = p±2- For m± ~ 3GeV, 5E = lOGeV, coshy ~ 100, the relative change in the pair 
spectrum is estimated to be F'(p±)/ F(p±) ~ 1 — 0.1(e~ y — 1). 

4. SUMMARY 

We calculated the lepton pair spectra from open charm decay in two different shad- 
owing scenarios. By scaling the ratios for different mass ranges according to eq.([5]), we 
showed that the dilepton rapidity dependence of those ratios on xa reproduces well the 
underlying gluon shadowing function defined in eq. ([]]). Finally we showed that the mea- 
surements required to extract the gluon shadowing are experimentally feasible at RHIC. 
The above analysis depends on the validity of our basic assumption regarding the hard 
fragmentation of charm quarks, and this can be checked explicitly via the single inclusive 
measurements of D production. 

We conclude by emphasizing the importance of determining gluon shadowing in p + A 
to fix theoretically the initial conditions in A + A. In A + A the open charm decay is 
regarded as an annoying large background that must be subtracted to uncover the thermal 
signal. In p + A that background becomes the signal needed to determine the incident 
gluon flux in A + A. The continuum charm dileptons in p + A at RHIC are likely to 
provide a unique source of information on the low xa nuclear gluon structure at least 
until HERA is capable of accelerating heavy nuclei. 
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